Off equilibrium response function in the one dimensional random field Ising model 
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I. INTRODUCTION 



A thorough numerical investigation of the slow dynamics in the d = 1 random field Ising model 
in the limit of an infinite ferromagnetic coupling is presented. Crossovers from the preasymptotic 
pure regime to the asymptotic Sinai regime are investigated for the average domain size, the au- 
tocorrelation function and staggered magnetization. By switching on an additional small random 
field at the time t m the linear off equilibrium response function is obtained, which displays as well 
the crossover from the nontrivial behavior of the d — 1 pure Ising model to the asymptotic behavior 
where it vanishes identically. 
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A phase ordering system, such as a ferromagnet quenched below the critical point, offers the simplest example of slow 
>^ [ relaxation with many of the interesting features observed also in more complex glassy systems jjj . At the core of this 
£^ ■ phenomenology is the separation of the time scales of fast and slow variables. In a domain forming system equilibrium 
is rapidly reached in the interior of domains while the interfacial degrees of freedom remain out of equilibrium for 
a time which diverges with the size of the system. In these conditions the order parameter autocorrelation function 
splits into the sum 

G(t,t w ) = G st (t-t w ) + G arS (t/t w ) (1) 

where G st (t — t w ) is the stationary time translation invariant (TTI) contribution due to fluctuations in the bulk of 
domains and G && (t/t w ) is the aging, or scaling contribution |^|,p| originating in the off equilibrium fluctuations. In 
the time scale over which G s t(t — t w ) decays to zero G &g {t/t w ) stays practically constant at M 2 , where M is the 
equilibrium value of the order parameter, while the decay of G ag (t/t w ) takes place for much larger time separations. 

[ Using G(t, t) = 1 for spin variables, this implies G st (0) = 1 — M 2 . 

A structure similar to ([!]) is observed also in the linear response at the time t to a random external field switched 

' on at the earlier time t w 

1 ■ X(t,t w ) =Xst(t ~t w )+ Xag(*,*«) (2) 

t3 : 

where the stationary contributions Xst(t — t w ) and G s t(t — t w ) are related by the equilibrium fluctuation dissipation 
O ■ theorem (FDT) 



T X st (t-t w ) = G st (0) - G st (t-t w ) (3) 

and X&g(t>~tw) is the off equilibrium extra response. Notice that (^) may be read as the statement that Xst{t — t w ) 
depends on time through G s t(t — t w ). Mean field theory for glassy systems predicts jl[ that in the asymptotic time 
region this holds also off equilibrium with x(t, t m ) = x(G(t, t w )). If this is the case and if lim^oo %(t, t w ) = x cq then 
static and dynamic properties are connected || by the relation 

-T^m) -*m (4) 



dG 2 



G=q 



where ^oq arid P e q(?) are the linear susceptibility and the overlap probability distribution |pj in the equilibrium state. 
Applying the scheme to phase ordering we should find 
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since the overlap function in the equilibrium state is given by the delta function H. On the other hand ([|) may be 
rewritten as T\st{t — t w ) = [1 — (G s t(t — t w ) + M 2 )] and, if t w is sufhciently large, parameterizing time through the 
full autocorrelation function we have 

^ ,~ J 1 - G , for M 2 < G < 1 . . 

1 Xst{G)-^ 1 _ M2 ,forG<M 2 (6) 

which gives — T^^Xst(G) ) = 5{q — M 2 ). Hence, for (Q) to hold one must necessarily have 



G=q 



lim X (t,t w )=Xst(G) (7) 



and lim^^oo x&g(t, t w ) — 0. Numerical simulations for the Ising model in d = 2 and d — 3 as well as analytical 
results for the spherical model || indeed show evidence that (0) is asymptotically satisfied. More precisely, one expects 
P,[To[ the off equilibrium contribution to the response function to scale as Xag(t, t w ) — t~ a x ag (t/t w ) with a — 1/2 on 
the basis of the argument that Xag(£, t w ) is proportional to the interface density pi(t) ~ where Lit) ~ t 1 / 2 is 

the typical domain size and z — 2 for non conserved order parameter dynamics || . 

Motivated by analytical results [[fl| for the one dimensional Ising model which, instead, give a — 0, recently a 



detailed study of the behavior of the response function has been undertaken 
the exponent a on space dimensionality. This is best explained in terms of t 
single interface and defined by 



finding an interesting dependence of 
re effective response Xes(t,t w ) due to a 



Xag(*, tw) = pl(t)XeK(t, tw)- (8) 

In this form the behavior of X^g(t,t w ) appears to be determined by the balance between the interface loss due to 
coarsening and the response associated to a single interface. If one requires Xag(Mto) ~ Pi(t) clearly Xeff(t,t w ) must 
be a constant. In Ref. |12] it was found that for an Ising system this is the case only for d > 3, while for d < 3 there 
is the power law growth 

Xeff(t,i w )~ (t-t w ) a (9) 

with numerical values for the exponent at different dimensionalities compatible with a = (3 — d)/4. At d — 3 the 
power law is replaced by logarithmic growth. It has been conjectured that this dimensionality dependence of a is 
the outcome of the competition between the curvature of interfaces and the external perturbing field in the drive of 
interface motion. According to this picture d = 3 is the dimensionality above which interface motion is dominated by 
the curvature, while below d = 3 the external field competes effectively with the curvature. The more so the lower is 
the dimensionality, until in d = 1 interfaces reduce to pointlike objects driven only by the external field. When this 
happens, the rate of growth of the single interface response matches exactly the rate of loss of the interface density, 
with a = 1/2 and a = 0. Then, Xa g (A t w ) does not vanish as t w — > oo and the second term in (||) contributes to x(G) 
producing the violation of (^). In short, dimensionality acts as the control parameter which allows to modulate the 
competition between the two opposing mechanisms driving interface motion and summarized by 

/ \ ,ford>3 . . 

fl= (%! ,ford<3. (10) 

In this paper we investigate a generalization of this phenomenon occurring in the framework of the random field 
Ising model (RFIM), which further clarifies how the overall behavior of Xag(ijiju) originates in the interplay between 
the rate of growth of the single interface response function and the rate of loss of defect density. Here we specialize to 
the d — 1 RFIM. In this case domain walls perform random walks in a random potential of the Sinai type Jl3| . The 
crucial feature of this walk, from our perspective, is the crossover between the preasymptotic regime, characterized by 
free diffusion as in the pure system, and the asymptotic regime characterized by Sinai diffusion. Measuring Xes(t,t w ) 
and pi(t), in the preasymptotic regime we shall find 

Xefl{t,t w ) >pj\t) (11) 

as in the pure system. Namely, as long as diffusion takes place in a flat landscape the rate of growth of Xeff(^i t w ) does 
not go below the rate of loss of interfaces. However, going over to the asymptotic regime the landscape becomes rugged 
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and activated processes do play a role. When this happens Xes(t,t w ) slows down with respect to pi(t) eventually 
reaching 

Xeff(*,0 ~ Pj^it) (12) 

which leads to the vanishing of Xag(t, t w ). Through this mechanism the validity of (||) is restored in the d = 1 RFIM. 



II. UNPERTURBED SYSTEM 



In the following we consider the d = 1 Ising model with Hamiltonian 

N N 

H[<Ji] = -jy~] (TjC7 i+1 - ^ h i a i ( 13 ) 
i=l i=l 

where J > is the ferromagnetic coupling and hi — ±ho is an uncorrelated random field with expectations 

E h (h i h j ) = h 2 S ij . li4j 

Here hi is not a perturbing field, so ho does not have to be small. We will consider the two different dynamical 
evolutions taking place when the system is quenched to the final temperature T starting from the two initial states: 

i) the system is prepared into a spin configuration containing a single interface 

<7i = sgn(?) (15) 

ii) the system is in equilibrium at infinite temperature with the uniform measure 

VoWi] = 2~ N . (16) 

In order to have a phase ordering process the equilibrium state at the final temperature of the quench must be ordered. 
In the d = 1 RFIM this is not possible even at T = 0, since the size of the ordered domains is limited by the Imry Ma 
length Ljm = So, we take the limit |l2]| of an infinite ferromagnetic coupling J — > oo in order to have Ljm = co- 
in this case the equilibrium state is the mixture of the two ordered states 

7%h = v^[-j = - + + !) 

i i 

for any finite temperature T. Notice that in the corresponding pure states M 2 = 1 which implies x cq = and 
P C q(q) — S(q — 1). Furthermore, with J — oo thermal fluctuations within domains are suppressed and dynamics 
is reduced to interface diffusion. With the initial condition (|l5|) we have the diffusion of the single interface and 
with the initial condition (^|) the diffusion-annihilation process of the set of interfaces seeded by the initial condition. 
Characteristic lengths of the process are in the first case the root mean square displacement (RMSD) R(t) of the single 
interface and in the second case the average distance between interfaces (or average domain size) L(t). Furthermore, 
since the typical size of the potential barrier that a walker must overcome after traveling a distance I is I 1 / 2 ho, there 
remains defined another characteristic length 



L > - UJ (18) 

as the distance over which potential barriers are of the order of magnitude of thermal energy This length is important 
because separates the distances smaller than L g , over which diffusion takes place as in the pure system with R(t) or 
L(t) ~ i 1 / 2 , from the distances larger than L g , over which diffusion is dominated by the random potential and it is of 
the Sinai type |Q with R(t) or L(t) ~ (logi) 2 . Clearly, the limit of the pure system corresponds to L g — oo. 

Considering different temperatures T in the presence of random fields of different strength h®, we have let the 
system evolve with initial conditions ([l5]) and (E6|). One time unit is defined equal to one flip attempt per spin, on 
the average (spins are updated in random orderJT^ We use Metropolis transition rates, that is the probability to flip 
the spin is pmp = min(l, exp(— yA7i)). In the case of initial condition (p"5|), we have found that the RMSD of the 
single interface depends on T and ho through L g satisfying (Fig. Ill) the scaling relation 
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R{t,L g ) = L g K[ — ) (19) 



with 



W~jf ,2 f rX<< ! (20) 

v ; [ (logx) 2 for x > 1 v ; 

which displays the crossover from the pure regime to the Sinai regime. A completely analogous behavior is obtained 
in the second case for the typical domain size L(t), measured as the inverse density of interfaces pj 1 ^). We find 
(Fig-!) 

L(t,L g )=L g c(j^j (21) 

where the scaling function C{x) obeys the same limiting behaviors as in (p0[). 

Next, let us consider the autocorrelation function G(t, t w , L g ) — jj-E^ {^2~[o'i(t)ai(t w ))h} where the angular brackets 
denote the average over thermal noise for a given realization of the field [hi] . Due to the absence of thermal fluctuations 
within domains G s t{t—t w ) = and the autocorrelation function is entirely given by the aging component which satisfies 
(Fig. ||) the scaling relation 



( L(t) 

\L(t w y 



G ae (t,t w ,L g ) = F(- r ^,z) (22) 



where z = M^il. The shape of the scaling function F(x, z) is known in the limits of the pure system (z = 0) and of 

the Sinai regime (z = oo). In the first case from ( ^l|) follows x = \pf- and from the exact solution |l5) of the Glauber 
dynamics for the Ising chain 



F(x, z = 0) = - arcsin I - ~ ., ) . ( 2>i i 



The second case is obtained taking t w so large that In L g can be neglected with respect to In t w . This yields x — ( -^—^J 



and flic 



F{x,z = oo) = -^=-±-. (24) 

O-y X OX 

In the intermediate cases with finite values of z one expects 

( , _ J F(x, 0) ,forx-Kl/z . . 

1 ' J ~\F(i,oo) ,forx-l>l/z 

where the condition on x can be rewritten more transparently as L(t) — L(t w ) <C L g or L(t) — L(t w ) 3> L g . While 
our data in Fig. |^ show clearly that (]2^) is obeyed for x — 1 1/z, the check of the crossover to the Sinai regime 
requires simulation times too long for what we can achieve. In any case, the behavior of the data under variation of z 
in Fig. H shows clearly the onset of the crossover. Notice that (^5|) means that even if the shortest time t w is chosen 
inside the Sinai regime, namely after the potential has developed a rough landscape, for displacements up to L g in 
the bottom of the potential valleys interface diffusion takes place as in the pure system. 
Finally, let us consider the behavior of the staggered magnetization 



M(t,L g ) = J ¥o E h 



(26) 



Essentially, this quantity gives the behavior of (minus) the magnetic energy per spin. At t = the spin and field 
configurations are uncorrelated and M(t = 0, L g ) = 0. As time evolves one expects the spins to correlate with the 
field producing growth in M(t, L g ). Indeed, M(t, L g ) displays (Fig. ||) growth in the preasymptotic regime L(t) <§; L g , 
which however is followed by an intermediate regime with a constant plateau for L(t) ~ L g and then by a decrease 
toward zero for L(t) ^> L g . In the intermediate and asymptotic regimes one has 



4 



M{t,L g )=M(^Pj (27) 

with M.(x) ~ x~ x / 2 for x ^> 1. In order to understand how this comes about it is useful to look at M(t,L g ) as 
made up by separate contributions associated to the single interfaces whose number decreases as time goes on. This 
is formalized by writing 

M(t,L g )= Pl (t)M eS (t,L g ) (28) 

which defines M c g(t, L g ) as the effective staggered magnetization associated to a single interface. Using ( p7| ) then we 
have 

M eS (t, L g ) = LgMcH (^Tj (29) 
with (Fig. ||) the function A4 e g(x) obeying 

The interpretation of M e s(t, L g ) as the single interface contribution to the build up of magnetization can be substan- 
tiated by measuring M(t^L g ) in the process with the initial condition (p"5|). In this case, defining M s - ms i e (t,L g ) — 
NM(t,L g ), we find (Fig. g) 

M single (t, Lg) = L g M singlc (^j^J (31) 

where A^ s ingie(^) displays the behavior (|30|). From this we may draw the following conclusions: i) when in the system 
there is only one interface, the spin-field correlation grows with time. Therefore, this is not a local effect involving only 
the alignment of the pair of spins at the interface with the local field. Rather, it is a large scale effect involving the 
optimization of the interface position with respect to the entire field configuration, it) The magnetization growth takes 
place with different time laws in the preasymptotic and asymptotic regimes. Hi) The total magnetization behavior of 
Fig. ^, when multiple interfaces are present, then is just due to the fact that in the preasymptotic regime the rate 
of growth of the single interface magnetization balances the rate of loss of interfaces, while in the asymptotic regime 
interfaces do disappear faster than the growth of the single interface contribution to the magnetization. 



III. RESPONSE FUNCTION 



Let us now consider what happens if at the time t w > an additional random field ti = ±eo uncorrelated with hi 
and with expectations 



E e (ei) = n 
EeleiCj) = e 2 5 k 



(32) 



is switched on. We take £q <^ ho and we are interested in the linear response function with respect to the e-perturbation 
given by 



T 

Tx(t,tw,L g ) = —^E h E € 
iv e n 



J2(vi(t))h^ t 



(33) 



where the external field h acts from t = and undergoes the change from h to h + e at t w . 

Due to the suppression of thermal fluctuations within domains enforced through ,/ = oo, the stationary contribution 
Xst(t — tto) in (0) vanishes identically. Therefore, the above response function is entirely constituted by the aging term 
Xag(i, t w ). The interest, then, is focused on the scaling properties. As in the case of the autocorrelation function, we 
find (Fig. 0) that this quantity obeys the scaling form 



T X (t,t 



L(t) 

L(t w ) 



,z 



(34) 
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For z = the pure system response function pi[ is recovered 

x(cc, z = 0) = arctan \J x 1 — I (35) 

which is responsible of the violation of (f|) in the pure Ising chain p2|| . In fact, from the equilibrium state ( |l7| ) we 
ought to have P e q(q) = ^(<? — 1) which is consistent with (Q), as explained in the Introduction, if only Xst(G) given 
by ( p|) e nters into and the limit M 2 — ► 1 is taken after differentiation. Instead, eliminating x between ( |35| ) 

and (|23| ) one finds 

X(G, z = 0) = — arctan[\/2 cot(^G)] (36) 

which spoils (|J). 

With z > there is a crossover. For values of x up to x — 1 ~ 1/z Fig. [?] shows that x(x, z) behaves as in the pure 
case, on the basis of the same argument used for the autocorrelation function. For larger values of x the data show 
that x(x, z) levels off and then decreases. This is clearly displayed also in the plot (Fig. ||) against the autocorrelation 
function. The mechanism responsible of this behavior is the same discussed in the previous Section for the staggered 
magnetization. Let us look at the effective response of a single interface Xeff(i> t w ) defined in (||). We find 

T XeS {t, t w ,L g ) = L(t w )x c s (^ttt, z] (37) 



L{t w y 



with the scaling function displaying (Fig. y) the behavior 



where Xes(x, z = 0) = xx(x, z = 0). ^From this follows ( |IT| ) in the preasymptotic regime and (|l2|) in the asymptotic 
regime, which account for the crossover of the response function in Fig. hj in terms of the balance between the rate of 
growth of the single interface response and the rate of loss of interfaces. Hence, for z > eventually X (x, z) vanishes 
and in the limit z — > oo one expects 

X(x, z = oo) = 0. (39) 

Therefore, for any finite ho the validity of (Eh is restored asymptotically. 

Lastly, it remains to make sure that Xcfil x i z ) may be identified with the response function of a single interface. 
This we have done by measuring the response function Xsingic(i, t w , L g ) in the process with only one interface in the 
initial condition and finding (Fig. [Tc]) a behavior quite close to the one of Fig. g. 



IV. CONCLUSIONS 

One of the hypothesis for the validity of relating static and dynamic properties is that the large time limit of 
x(t, t w ) reaches the equilibrium value Xoq- It is then clear from (Q) that for this to be true in the phase ordering process 
Xag(^to) must vanish, since this is an intrinsically out of equilibrium contribution. In other words, the existence of 
the interfacial degrees of freedom which do not equilibrate must not play a role at the level of the response function. 
Indeed, X a.g(t,t w ) normally does to vanish. However, there are exceptions in special cases and the study of these 
is quite instructive since it allows to gain insight into the properties of Xa.g(t,t w ). One dimensional systems with a 
scalar order parameter and frozen bulk thermal fluctuations do make such a special case. In the d — 1 pure Ising 
model, due to the pointlike structure of domain walls, the minimization of magnetic energy introduces a bias in 
the diffusion of interfaces which leads to a nonvanishing Xng(t,t w ) through (fy]). Elsewhere we have analyzed 
how ([ll]) ceases to hold in going from d = 1 to d > 1, due to the extended nature of interfaces. In that case the 
minimization of magnetic energy is hindered by the competing need to minimize the curvature of interfaces and 
Xag(t,t w ) asymptotically disappears. In this paper we have investigated a different mechanism altering the delicate 
balance ( |IT| ) between the gain and loss of contributions to the response function, which does not require the passage 
from pointlike to extended defects. In the d — 1 RFIM it is the gradual roughening of the landscape that slows down 
the minimization of the magnetic energy with respect to the growth law of the domain size and which eventually 
yields an asymptotically vanishing Xag (£,£«>)• Therefore, the overall picture which comes out is that a coarsening 
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system is always out of equilibrium in the sense that there are always interfaces around, each one producing the out 
of equilibrium contribution Xeff{t,t w ) to the response function, independently from dimensionality or the presence of 
quenched disorder. These are elements which affect the finer details, like whether Xcs(t,t w ) is a constant or how fast 
it grows. Where these properties come into play is in putting together the contributions of all the separate interfaces 
through Xag(i, tw) — pi(t)xes(t, t w ) which, now depending on dimensionality or quenched disorder, may then produce 
either a vanishing or an asymptotically persistent net result. 

Finally, a comment should be made about the non commutativity of the order of the limits ho — > and t w — > oo. 
If the limit ho — > is taken first, the linear response function of the pure Ising model is obtained and, as discussed 
above, (Q) is violated. Instead, if tw — > oo is taken first one finds the asymptotic linear response function (^) of 
the RFIM which is consistent with (^). However, if the limit ho — > is taken next, the response function sticks to 
x(x, z — oo), which is not the linear response function of the pure model. 
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FIG. 1. Rescaled RMSD R(t)/L g versus rescaled time t/L 2 g , for a single interface and L g — 1,4,25. The fit parameters are 
a = 1.33, b = 0.274 and c = 1.33. 
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FIG. 2. Rescaled average length of domains L(t)/L g versus rescaled time t/L 2 g , for multiple interfaces and L g — 400, 1600. 
The fit parameters are a = 3.45, b = 0.271 and c = 15.8. 
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FIG. 3. Autocorrelation function G(t,t w ) versus the ratio {L(t) — L(t m ))/ L(t w ), for L g — 400, 1600, and z = 0.1, 1, 4, with 
z — L(t w )/L g . Solid lines are the exact results for z — > and z — > oo. 
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FIG. 4. Total staggered magnetization M(i) versus rescaled average length of domains L(t)/L g , for L g — 400, 1600 




FIG. 5. Rescaled effective staggered magnetization M e g(t)/L g versus rescaled average length of domains L(t)/L g , for L g 
100. The fit parameters are a = 0.68 and b — 1.7. 




FIG. 6. Rescaled single interface staggered magnetization M S i ng i e (£)/ L g versus rescaled RMSD R(t)/L g , for L g — 1,4,25. 
The fit parameters are a = 0.75 and b — 1.15. 
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FIG. 7. Total linear response Tx(t,t w ) versus the ratio (L(t) - L{t w ))/L{t w ), for L g = 400, 1600, and z = 0.1, 1, 9. The 
solid line is the exact result for z — > 0. 
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FIG. 8. Total linear response T\{t, t w ) versus autocorrelation G(t,t w ), for L g — 25,100 and z = 1,4,9,16. The solid line is 
the exact result for z — > 0. 
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FIG. 9. Rescaled effective linear response Tx e s{t, t w )/L(t w ) versus the ratio (L(t) — L(t w )) / L(t w ) , for L g = 25, 100, and 
z = 1, 16. The solid line is the exact result for z — ► 0. 



11 




10 2 10 1 1 10 10 2 



R(t,t w )/R(t w ) 



FIG. 10. Rescaled single interface linear response Tx S i ng i e (£, t w )/R(t w ) versus the ratio R(t,t w )/R(t w ), where R(t,t w ) and 
R(t w ) are respectively the RMSDs between times t w and t, and between times and t w , for L g = 4, 25, and z = 1, 9. The fit 
parameters are a — 0.75, b(l) = 1.1, 6(9) = 0.35. 
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